Scrambling eggs: Meiotic drive and the evolution of 
female recombination rates. 
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Abstract: Theories to explain the prevalence of sex and recombination have long been 
a central theme of evolutionary biology. Yet despite decades of attention dedicated to the 
evolution of sex and recombination, the widespread pattern of sex-differences in the recom- 
bination rate is not well understood and has received relatively little theoretical attention. 
Here, we argue that female meiotic drivers - alleles that increase in frequency by exploiting 
the asymmetric cell division of oogenesis - present a potent selective pressure favoring the 
modification of the female recombination rate. Because recombination plays a central role in 
shaping patterns of variation within and among dyads, modifiers of the female recombination 
rate can function as potent suppressors or enhancers of female meiotic drive. We show that 
when female recombination modifiers are unlinked to female drivers, recombination modi- 
fiers that suppress harmful female drive can spread. By contrast, a recombination modifier 
tightly linked to a driver can increase in frequency by enhancing female drive. Our results 
predict that rapidly evolving female recombination rates, particularly around centromeres, 
should be a common outcome of meiotic drive. We discuss how selection to modify the 
efficacy of meiotic drive may contribute to commonly observed patterns of sex- differences in 
recombination. 
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1 Background 



Recombination plays a critical role in both the production of viable gametes and population 
genetics processes. Structurally, chiasmata - the physical manifestations of crossovers - 
generate the tension between homologs often needed to ensure proper segregation during 
meiosis I. The structural role of chiasmata is likely the mechanistic underpinning of the 
requirement in most species of at least one recombination event per chromosome. In addition 
to the structural role of chiasmata, recombination also plays an important role in population 
genetic processes by generating novel haplotypes within populations. The production of 
novel haplotypes entails both the creation and separation of beneficial alleles (or allelic 
combinations). The balance between these o pposing outcomes shape the adaptive value o 



recombination, a topic of much research (e.g 



1996 



Otto and Lenormand 



2002 



Eshel and Feldman 



Keightley and Ott 



2006 



Re combination rates vary between species, i ndividuals, and sexes ( 



1988 



2005 



Burt et al 



1991 



Lenormand 



Coop and Przeworski 



2003 



1970 



Fe 



Barton 



BellI. 



Lenormand and Dutheil 



LDMAN et al 



2009L 



1982 



Trivers 



2005 



Lorch 



20071 ). We address the evolution of sex-difference in the 
recombination rate (i.e. heterochiasmy). We first characterize general patterns of sex- 
differences in recombination rates. After reviewing the major theories invoked to explain 
the evolution of heterochiasmy, we introduce the major argument of this manuscript - that 
the common and consistent sex-difference in the operation of meiotic drive can favor the 
evolution of heterochiasmy. We conclude by considering the implications of our model in 
light of our current understanding of the genetic basis of recombination modification. 



1.1 Observations 

Four broad patterns describe sex differences in recombination rates. We describe these 
patterns and their consistency below. 

The achiasmatic sex is heterogametic (The Haldane- Huxley Rule): When recombination 
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is absent in one sex (i.e. it is achiasmatic), that sex is ne arly always the heterogame t ic sex 
(i.e. the sex beari ng heteromorphic sex chromosomes 



Burt et al. 



Haldane. 


1922; 


Huxley. 


1928 



19911 ) . This obser vation represents m ore than 25 evolutionary independent 



19911 ). and is observed under both male (e.g. 



origins of sex-specific achiasmy (IBurt et all 
Drosophila) and female (e.g. Lepidoptera) heterogamy, with very few known exceptions 



( IDavies and Roderick 



2005|). 



Females recombine more than males: When both sexes recombine, t 



re female recombi- 



nation rate oi 


'ten exceeds the male rate ( 


Bell. 


1982; 


TRiVERS: 




1988 


Lenormand. 


2003; 


Lenormand and Dutheil 


2005; 


LORCH 


• 


2005 


)• 



Burt et al. 



1991 



tern in Figure [TJA. by plotting the sex difference in recombination rates (measured by auto- 
somal map lengths or chiasmata counts) in a number of taxa. We note that the grouping of 
our taxa is somewhat arbitrary; nonetheless, clear trends across groups in the degree of sex 
difference in recombination rate are apparent. 

The pattern of higher female recombination rates is quite broad, occurring in animal 
species with XY, ZW, and environmental sex-determination. However, there are many ex- 
ceptions (e.g. marsupials, some grasshoppers, and newts) suggesting that the process of 
recombination is not mechanistically constrained towards high er rates in females and that 
the ra tio of male to female rates can evolve quite rapidly (see ILenormand and Dutheil 



( 120051 )). In plants, there is no clear trend towards higher recombination rates in female 
meiosis (Figure [TJ\); however, when outbreeding angiosperms are considered separately, re- 
combination rates are on average slightly higher in female meiosis than in male meiosis 



( ILenormand and DuthehJ . 



20051 ) 



Females recombine at relatively higher rates near centromeres: After controlling for the 
genome-wide sex difference in recombination rate, females recombine more often near cen- 
tromeres than do males, while males recombine re l atively more near telomeres. This patter n 



has been noted in fish (Sakamoto et al 



humans ( IBroman et al. 



1998 



2000 



Kong et al 



2002 



Singer et al. 



Clark et al 



2002 



Reid et al 



200711 



2010), dogs ( IWONG et al 
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A) Comparative pattern 



B) Genomic pattern 




Figure 1: Genome wide and regional sex-differences in recombination rates. A) 

The difference between female and male recombination rates ( o from linkage maps, excluding 
known sex-chromosomes, x from chiasmata counts) divided through by the sex-averaged 
rates. Points above the dashed line indicate higher rates of recombination in females than in 
males. * indicates p < 0.05 using a two-tailed sign test, without correcting for multiple tests 
or phylogeny, and ignoring ties. B). Sex-standardized recombination rates across the human 
genome. The sex-standardized rate equals the local recombination rate in a given sex (male 
and feamlc), divided by the average recombination rate in that sex. The x-axis indicates 
the position of the focal genomic region (.2% of a chromosome arm), divided by the length 
of the chromosome arm. Data are presented from all metacentric human autosomes. Lines 
represent a lowess smoothing of these points. 
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20 id ), and mice (IShifman et al. 



tions (e.g. opossums 



2006 



Samollow et al. 



Paigen et al. 



2004 



20081 ) , although there are some excep- 



20071 ). Ut ilizing data 



rom a recent fine-scale 



2010l ). we display an 



analysis of sex-specific recombination rate in humans flKoNG et al. 
example of this pattern in Figure IB. 

This pattern is not an obvious consequence of different requirements for segregation in 
males and females, as it is thought that recombination near centromeres does not contribute 
to cohesion of homologous chromosomes. Furthermore, recombination events too close to 
the centromere can generate problems during meiosis, such as an increase in the rate of 
precocious separation of sister chromatids, potentiall y leading to aneuploid gametes (i.e. 



gametes with an abnormal number of chromosomes - 



Rockmill et al. 



20061). The higher 



rates of recombina tion close to centromere in aneuploid transmissions ( IMay et al. 



1990 



Lamb et al. 



20051 ) suggests that recombination close to the centromere may actually incur 



fitness costs due to errors during meiosis. 

Genetic control of the recombination rate is often sex specific: The process of gametoge- 
nesis is fundamentally different in males and femal es. Consequently, the genetic control of 



20051 ). Perhaps 



many aspects of meiosis differs between the sexes (IMorelli and Cohen! 
due to sex differences in meiosis, alleles that influence the recombination rate in one sex often 
have no influence (e.g. the 17q21.31 inversion region in humans), or the opposite effect (e.g . 



RNF21 2 in humans) on the r e combi nation rate in the other sex (IKong et al. 



2004 . 



2008|). 



In fact, IFledel-Alon et all ( 1201 ll ) show that there is no detectable heritable intersexual 
correlation in the recombination rate despite additive genetic variance in both sexes. 

Although mechanisms governing sex differences in the recombination rate are not well- 
characterized, one candidate is the length of the synaptonemal complex (SC), a struc- 
ture composed of cohesins and other proteins involved in meiosis, which stabilizes connec- 
tions between homologs and facilitates crossing over. SC length is p ositively associatec 



with the recombinat i on ra t e, and is longer in fema 



Tease and Hulten 



2001 



Dumont and Payseur 



es than in males ( Lynn et al. 



2002 



20111 ) . Regardless of mechanism, the 
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consequence of this pattern is clear - there is ample opportunity for sex-specific recombination 
rates in one sex to evolve independently of those of the other sex. 



1.2 Current theories 



Des pite the large body of work on t 



sec 



Otto and Lenormand 



i e evolution of se x and recombination (for recent reviews 



2002 



Barton 



20091 ). the evolutionary forces responsible for 



the obser yed patterns of 



(but see 



Trivers 



h eterochiasmy have r e ceived curious ly little theoretical attention 



1988 



Lenormand . 



2003 



LORCH 



20051 ). Below, we summarize the 



current of the evolution of heterochiasmy and review their plausibility. 

We note that these theories are not mutually exclusive, and it is unlikely that any one 
theory will explain all the observations presented above. Further, our model (like many of the 
other theories of the evolution of heterochiasmy) does not argue that selection directly favors 
the evolution of male and female rates in opposite directions. Rather, in our model, selection 
favors a modification in the female rate, with little or no direct selection on the male rate. In 
both our model and others, heterochiasmy results from the sex-specificity of recombination 
modifiers and/or additional selective constraints (such as stabilizing selection on the sex- 
averaged rate) which are not made explicit. For example, if too much recombination is 
generally deleterious due to the costs of ectopic recombination, but circumstances favor an 
increase in female recombination rates, modifiers that increase female rates but not male 
rates may be favored, leading to the evolution of heterochiasmy. 

Sex chro mosome pl e iotrop y: B ased on their ob servations that the achiasmatic sex is het- 



erogametic, 



HaldaneI (119221 ) and 



Huxleyi (119281 ) proposed that heterochiasmy may evolve 



as a pleiotropic consequen ce of select ion to suppress recombination between heteromorphic 



sex chromosomes (see also 



Nei 



1969)- Although suppression of recombination between sex 
chromosomes could well explain the Haldane-Huxley rule, it cannot fully explain the quan- 
titative variation in autosomal recombination rates between the sexes, as they span modes 
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of sex determination, nor can this model explain variation across chromosomal regions. 

Maternal aging: Physical connections between homologous chromosomes are necessary 
for proper segregation during meiosis. In many species, these physical connections are formed 
by chiasmata and therefore, one recombination event per chromosome is generally required 
for proper meiosis. In species where female meio sis is arrested, phy sical connections be- 



Lamb et al. 



20051 ). and thus additional 



tween chromosomes may degrade with time (e.g. 
chiasmata may function to stabilize chromosomes across the metaphase plate. 

According to the maternal aging theory, elevated female recombination rates provide more 
physical connections between chromosomes, ensuring that oocytes will segregate correctly 
after years of insults accumulated during meiotic pachytene arrest. This theor y is supported 



by th e finding that women with higher recombination rates have higher fertility ( IKong et al. 



2004 ) and that the viable gametes of older mothers have a higher number of crossovers, 



consistent with the idea that selection aga inst those gametes with too few crossove rs increases 



throughout a mother's reproductive life (IKong et al 



2004 



Coop et al 



20081 ). However, 



this theory cannot easily explain heterochiasmy in organisms in which females create eggs 
throughout their lives. Additionally, this theory cannot easily explain the spatial pattern of 
recombination in females, since elevated female rates are often accomplished by chiasmata 



placed in locations t 

Sexual selection: 



lought not to st abilize chromosomes. 



Trivers 



(119881 ) argued that selection to preserve high fitness genotypes 



favors r ecombination suppr ession in t he sex with greater variance in fi tness. Both current 



theory (ILenormand 



20031 ) and data ( IBurt et al 



1991 



Mank 



20091 ) suggest that sexual 



selectio n cannot expla i n the evolution of heterochiasmy. Using multi-locus population genetic 



theory, 



Lenormand 



(120031 ) showed that sex- differe nces in select i on on dipl oid genotypes 



Burt et al. 


(1991 


) and 


Mank 


(2009) 



cannot generally favor the evolution of heterochiasmy. 
showed that the degree of heterochiasmy decreases with sexual-size dimorphism (a proxy for 
the strength of sexual selection) - an observation counter to predictions of the sexual selection 
hypothesis. 
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Haploid and psuedo-haploid selection: 



LenormandI (120031 ) showed that simple sex differ- 



ences in selection during the haploid life stage can favor the evolution of heterochiasmy. Like 
the sexual selection theory, the haploid selection model argues that the sex which produces 
the gamete (or gametophyte) with higher variance in fitness will recombine l ess. As sperm 



Lenormand 



Lenormand 



and p ollen can experience intense competition for fertilization opportunities 
( 2003 ) argued that males should in general have lower recombination rates. 
(120031 ) also proposed that selection on haploid components of diploid genotypes (e.g. selec- 
tion on epistasis in cis, or imprinted loci) can also favor the evolution of heterochiasmy; we 
call this the pseudo-haploid selection model. Although theoretically plausible, it is unclear 



whether the small numbers of imprinted 



sperm (e.g. 



Joseph and Kirkpatrick 



genes 



(Imprison et al. 



20051 ) or genes expressed in 



20041 ) are sufficient for the (pseudo-)haploid selec- 



tion theory to explain heterochiasmy in animals. Furthermore, a comparative study found 
no association between heterochiasm y rates and the inferred strength of sperm competition 



in eutherian mammals (IMank . 



20091 ). However, the absence of a haploid stage in the female 
gametes of animals and occasional haploid expression in sperm makes the haploid selection 
theory viable. 

The situation is more complex in plants, as due to the alternation of generat ions, there is 
haploi d expression in both the products of male and female meiosis. Nonetheless 



Lenormand and Duthe 



( 120051 ) argue that as in the majority of outcrossing plant species, selection is likely strong on 
male haploid products (due to pollen grain competition), the haploid-selection model could 
explain the female biased recombination rates in outcrossers. 

Meiotic drive: Below, we articulate the meiotic drive hypothesis. According to this 
model, sex differences in the operation of gametic drive presents a sex-specific selective 
pressure on linked and unlinked modifiers of the recombination rate. Other authors have 
proposed that sex differences in meiotic drive may of f er an opport unity for the evolution of 



sex-specific recombination rates ( ILenormand 
relation to these models below. 



2003; 



Haig 



20101 ) . We discuss our work in 



9 



1.3 Meiotic drive and the evolution of heterochiasmy 



Meiosis provides an opportunity for alternative alleles to compete for representation in the 
functional gametes of heterozygotes. Alleles that distort meiosis and gametogenesis in their 
favor (i.e. gametic drivers) often do so at the expense of individual viability or fertility. 
Therefore, although driving alleles can benefit by distort ing mei osis, individual selection 



generally fa vors Mendelian segregation (IEshelI e.g 



Haig 



3 e.g. 



1985 



but see 



Ubeda and Haig 



2005 



20101 ). creatin g a conflict be tween drivers and unlinked loci in the genome (i.e. 'the 



parliament of genes,' 



Leigh 



197lh . 



Gametic drivers exploit the system of Mendelian segregation by providing a transmis- 
sion advantage to their chromosome. Higher recombination rates make the ultimate chro- 



mosomal context of an alle 



Thomson and Feldman 



e uncertain, which can prevent the evolu t ion of drive systems 



1974 ; 



Charlesworth and Hartl 



1978 



Haig and Grafen 



19911 ). It is therefore thought that modification of the recombination rate can evolve as 
a mechanism to alter the efficiency of gametic drivers. Conceptually, this model holds for 
both male and female drive systems. However, gametic drivers are often sex-limited and dis - 



play sex differences in the mechanisms by which they operate (see 



Ubeda and Haig 



2005|). 



Currently, the implications of sex differe nces in meiot ic drive for the evolution of sex-specific 



recombination rates are unclear (but see 



Haig 



2010 



for an hypothesis) 



Male gametic drivers (e.g. Segregation Distorter in Drosophila and the t-haplotype in 
mice) u sually operate after meiosi s and are characterized by a two- locus damage-insensitive 



system (IWu and Hammer 



19911 ). When these loci are tightly linked, a damage-insensitive 



haplotype can incre a se in frequency, even if this haplotyp e decreases individual fitness (e.g. 



Prqut et al 



1973 



Charlesworth and Hartl 



19781 ). If the drive system imposes a 



cost, un linked recombination enhancers c a n be favored for the ir abili ty to disrupt the drive 



system ([Thomson and Feldman 



1974 



Haig and Grafen 



19911 ). However, modifiers 



of the recombination rate linked to and in phase with driving alleles can increase in fre- 
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quency if they decrease t 



l e rec o mbination rate between component s of a drive system 



(IThomson and Feldman 



1974 



Charlesworth and Hartl 



19781 ) . This later idea is 



supported by the observat ion that gametic drivers in males are often located in inversions 



fjWu and Hammer! 



19911 ) . which act to suppress recombination locally. 



Female meiosis creates a single haploid product, provid 
compete for representation in the egg (true meiotic drive 


mg an opportunity for al 


c 


es to 


Sandler and I" 


"lOVITSKI 


1957; 


Zwick et al. 


1999: 


Pardo-Manuel De Villena and Sapienza. 


2001b 


a ) . Recombina- 



tion plays a fundamental role in female drive systems, as it determines at which stage of 
meiosis alleles can compete with each other. Since non-sister centromeres segregate at 
meiosis one (hereafter MI), an allele that biases the outcome of MI in favor of its cen- 
tromere becomes over-represented in oocytes, so long as there is no (or an even number of) 
recombination (events) between driver and centromere (see Figure [2A.1 and 2A.2). The 



best-characterized cases of MI female drive are a subset o 



mammals (IPardo- Manuel De Villena and Sa 



in Mimulus guttatus 
and Henikoff (e.g. 



Fishman and Willis 



Malik et al. 



2002; 



2005 



ienza 



Robertsonian translocations in 



2001dT). and a cen t romer ic allele 



Fishman and Saunders 



Malik and Henikoff 



2009 



Malik 



2008h. Malik 



2009|) have 



argued for a broad role, throughout eukaryotic evolution, of female centromeric drive at MI 
in driving the rapid evolution of centromeric sequence and the proteins that bind them. 

Centromeres cannot drive during meiosis II (hereafter Mil), as they are paired with their 
sister chromosomes. However, with a single recombination (or odd number of) event (s) 
between a focal locus and the centromere, the products of the first meiotic dividsion (dyads) 



will be variable at that locus, providing an opportunity for Mil drive 



see Fi gure W) ■ 



Exa mples of Mil drivers include the Om locus in mice (IWu et al. 



maize (IRhqades and Dempsey 



2005|) and Ab-10 in 



19661 ). The structure of the Ab-10 haplotype (an inver- 
sion spanning multiple loci) highlights the importance of recombination in the evolution of 
female drivers. A combination of the alleles i n the Ab-10 system allows the c hromosome to 



drive during Mil (IRhqades and Dempsey 



1966 



Dawe and Cande 



19961 ). while a dis- 
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tinct allele at a nother locus in this compl ex alters the recombination rate b etween itself anc 



its centromere (IHiatt and Dawe 



1999|). 



20031 ). maximizing its ability to drive (IBuckler et al. 



In this manuscript we show that sex differences in meiotic drive can favor the evolution 
of heterochiasmy. With female meiotic drive, female recombination modifiers (unlinked to 
drivers) are favored when they decrease the efficacy of female drive. With segregating MI 
drivers, this corresponds to a female recombination enhancer, while female recombination 
suppressors decrease the efficacy of Mil drivers. By contrast, when recombination modifiers 
and drivers are in one tightly linked haplotype, female recombination modifiers that increase 
the efficacy of drive are favored (Table 1). 



1.4 Relation to previous models 

Numerous authors have explored cases in which recombination m odifiers are favored for their 



ability to break up systems of transmission r atio distortion (e.g. 



1974 



Charlesworth and Hartl 



1978|). 



'homson and Feldman 



Haig and Grafeni (jl99ll ) pointed out that in 



addition to breaking apart drive systems, recombination can decrease the efficiency of drive 
by making the identity of an allele's partner in a dyad uncertain. 

The po tential role of meioti c driv e in the evolution of heteroc h iasmy has received less 



attention. 



LenormandI (120031 ) and 



Lenormand and DutheilI (120051 ) briefly discussed 



male gametic drive systems, as a special case of the haploid selection model that may favor 
the evolution of heterochiasmy. As haploid chromosomes drive, the female meiotic drive 
model could also be considered as a special instance of the haploid selection model. However, 
we present our model as a distinct hypothesis because of its focus on recombination as a 
mechanism t o modi f y the efficacy of meiotic drivers. 



Recently, 



HaigI (120101 ) found that unlinked modifiers of the female recombination rate 



increase in frequency when they enhance the efficacy of drivers (the opposite of our finding, 
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below). In both our model and Haig's, drivers decrease individual fitness. However, under 
Haig's model, the cost is reflected in the genetic identity of the products of meiosis, in which 
the fertility of a female heterozygous for a driver is determined by the genetic composition 
of her dyads. Haig's model equilibrates when the cost of drive on fertility is balanced by 
the degree of distortion. At this equilibrium, a modification of the recombination rate which 
increases female fertility also results in an increase in frequency of the driver, leading to the 
counterintuitive result that recombination modifiers unlinked to driver benefit by increasing 
the efficacy of drive. In other words, in Haig's model, females produce a higher proportion 
of viable gametes by creating chromosomal configurations that lessen the ill effects of drive. 
Somewhat paradoxically, by doing so, unlinked modifiers act to increase the frequency of 
drivers, despite the long term population-level fitness consequences of the spread of drivers. 

In the model described below, we come to the opposite conclusion. The discrepancy 
between our models arises from different conceptions of how drive influences fitness. In our 
model an individual's viability depends on it genotype at the drive locus, but the viability 
of all dyads are equivalent. Therefore, the recombination modifier increases in frequency 
by decreasing the strength of drive. By contrast, in Haig's model, recombination directly 
influences female fertility by influencing the viability of eggs. Which model is more relevant 
to the evolution of heterochiasmy will depend on the mechanistic underpinning of the the 
cost of female drive. 

2 The models 

Since neutral and beneficial drivers will quickly sweep through a population, and will not 
present a genetic conflict, we focus on the evolution of drivers that decrease organismal 
fitness. We explore two models of drive. In the first model, the drive system consists of 
a single-locus that drives at either MI or MIL In our second model, drive is achieved by a 
two-locus system in which the ability of an MI centromeric driver to distort meiosis depends 
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Case 


When 


Linkage 


Prediction 


Reason 


1 


MI 


Unlinked, or 
repulsion phase 


r 9 > r d 


Female recombination enhancers dis- 
courage drive and hitchhike with high 
fitness, non-driving haplotypes. 


2 


MI 


Coupling 


r 9 < r<s 


Female recombination suppressors 
hitchhike with driving haplotypes. 


3 


Mil 


Unlinked, or 
repulsion phase 


r 9 < r«5 


Female recombination suppressors dis- 
courage drive and hitchhike with high 
fitness, non-driving haplotypes. 


4 


Mil 


Coupling phase 


r 9 > r 6 


Female recombination enhancers hitch- 
hike with driving haplotypes. 



Table 1: Relative male (r^) to female (r 9 ) recombination rates as predicted by the meiotic 
drive theory. 



on the genotype at a partially linked drive enhancer locus. 

For the single-locus drive model, we only care whether sister alleles are or are not sep- 
arated during the first meiotic division (i.e. if there there is an odd (1, 3, 5, • • • ) or even 
(0, 2, 4, • • • ) number of recombination event occur between the centromere and the drive 
locus) in drive heterozygotes. In our model, sisters are separated at the first meiotic division 
with probability r, or remain united with probability, 1 — r. 

With an even number of recombination events between drive and centromeric loci, there 
will be variation within, but not among dyads. This genetic difference among dyads presents 
an opportunity for drive during MI, but prevents drive during MIL By contrast, an odd 
number of recombination events between drive and centromeric loci produces genetically 
homogenous dyads, precluding the possibility of drive during MI, but facilitating Mil drive. 
Given an opportunity for drive, drivers are represented in ayii > 1/2 of the gametes produced 
by drive heterozygotes. Without an opportunity for drive, meiosis is fair. 
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We then introduce another layer of biological complexity - the dependence of MI drive on 
the two-locus genotype at centromeric and drive enhancer loci. In this two-locus system, the 
genotype at the drive enhancer locus influences the ability of the driving centromeric allele 
to distort segregation at MI. Without a drive enhancer, meiosis is fair. By contrast, the 
driving centromere is present in a,\ of gametes from double heterozygotes, and a<i of drive 
enhancer homozygotes heterozygous at the centromeric drive locus. As Mil drive is unlikely 
to depend on the genetic identity of the centromere, we do not explore the case of two-locus 
drive during meiosis II. 

Under this two-locus model, the outcome of multiple recombination events is somewhat 
complicated. Therefore, we make the restrictive assumption that only zero or one recombi- 
nation events occur on a tetrad between the centromere and the drive locus (i.e. complete 
crossover interference). Thus, in our two- locus model r is the probability of a single recom- 
bination event between the centromere locus and the drive locus, 1 — r is the probability of 
no recombination, and we assume no double crossovers. 

For each model, we contrast the case of sex-limited drive to results from a model in which 
both sexes drive. In these models, we introduce an allele that modifies the recombination 
rate between driver and centromere, without otherwise influencing organismal fitness. For 
models of single-locus drive, we compare the evolution of recombination modifiers in tight 
linkage with drivers to the evolution of recombination modifiers unlinked to drivers. 

3 Single-locus drive 

Consider a biallelic locus at which the driving allele, D, occurs in frequency, fjj, and the 
alternative allele occurs in frequency fa = 1 — fr>- As we assume that drive has a pleiotropic 
cost on individual fitness, heterozygotes and drive homozygotes have fitnesses wod ^ 1 and 
wdd < 1; respectively. Although we assume that this cost in suffered equally by both sexes, 
allowing for sex-specific costs does not change the qualitative picture. 
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|®@©(D 



Oppotunity for drive? 
Odd # recs Even # recs 



First Division (Ml Drive) No 


Yes 


Second Division (Mil Drive) Yes 


No 



Figure 2: Recombination during oogenesis and the opportunity for meiotic drive. 

Only one of the four products of female meiosis is included in the egg. Recombination is 
a critical determinant for the opportunity for drive because it partitions variation within 
and among the products of the first meiotic division (dyads). With recombination between 
the marker and centromere, there is variation within but not among dyads, presenting an 
opportunity for drive during Mil but not MI. When recombination occurs after the marker, 
there is variation among but not within dyads, presenting an opportunity for drive during 
MI but not MIL 
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The effective strength of drive acting in the i th meiotic division, Xya, is a function of 
both the recombination rate, r, and the strength of distortion, aui- With drive in both 
sexes, Xmi = «mi(1 — r ) + r /2, and Xmii = raun + (1 — r)/2, for MI and Mil drivers, 
respectively. With female-limited drive, the strength of MI and Mil drive in oogenesis is 
X M iq = a MI 9(l — rg) + rg/2, and X M ug = TOmiip + (1 — r 9)/2, respectively. Note that 
the effective strength of MI drive decreases with the recombination rate, while the effective 
strength of Mil drive increases with the recombination rate. Note further that the effective 
strength of female drive is independent of the male recombination rate. 

When drive has equivalent fitness effects in males and females, the mean population 
fitness, W, equals 

W = fdd + fDdWDd + fDDWDD (1) 

where /. denotes genotypic frequencies in newborns. If drive operates in both sexes, the 
frequency of a driver after selection and drive is equivalent in sperm and eggs, and equals 

JdMi ~ w [ ' 

When drive is female-limited, the recursive equation describing the frequency of drivers in 
sperm after selection (f' D(5 ) is 

,/ _ w DD f DD + w Dd f Dd /2 

Jd6 w { ■ ) 

and the recursive equation describing the frequency of drivers in eggs {f' D g) is 

w WDDfDD + ^ Ddf DdXy^g . . 

J D9M% ~ W 1 ' 

where genotypic frequencies equal f DD = f D QfD6, fod = Id? fdd + fdvfod, and f dd = f d gfdd, 
and the subscripts, 9 and 6, represent allele frequencies in female and male gametes, re- 
spectively (i.e. egg and sperm). With drive in both sexes, allele frequencies are identical in 
sperm and eggs, and therefore genotypes are in Hardy- Weinberg Equilibrium; however, with 
female-limited drive, Hardy- Weinberg assumptions are inappropriate. 
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In our appendix, we derive equilibrium frequencies of drive alleles in the absence of 
recombination modifiers by solving for fjj when Afjj = 0. Since a protected polymorphism 
at the drive locus facilitates the evolution of recombination modifiers, the existence of these 
equilibria is important for the models below. 

3.1 Single-locus drive - An unlinked recombination modifier. 

We now investigate the coevolution of alleles at the drive locus and alleles at a locus which 
influences the recombination rate. At the recombination modifier locus, alleles M and m 
occur in frequencies fu and f m = 1 — /at, respectively. The M allele additively alters the rate 
of recombination between driver and the centromere, has no direct effect on individual fitness, 
and is unlinked to the drive locus. The recombination rate between drive and centromeric 
loci is r, r + 8r, and r + 25r in mm, Mm, and MM individuals, respectively. In the case of 
drive in both sexes, we allow M to have equivalent effects in males and females. With single- 
locus, female-limited drive we only consider the influence of M on the female recombination 
rate, since our results are unaffected by r<j. We note that since M and D are unlinked, M 
does not influence the rate of recombination between itself and the drive locus. 

The population consists of four haplotypes (md, Md, mD, and MD). Equations describ- 
ing haplotype frequencies after selection, recombination, and drive, as well as many other 
derivations are presented in our appendix. 

To find the frequency of the recombination modifier, M, after one generation, we sum 
the frequencies of both haplotypes containing the M allele. With drive in both sexes, the 
change in frequency of M after selection, recombination, and drive equals 



where LD is the linkage disequilibrium between driving and recombination modifier alleles 
(i.e. LD = Jmd — /m/d)j and is measured in the gametes which united at random to form 
this generation. 




(3.a) 
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Equation (13. ap shows that a genetic association between drive and modifier loci {LD 0) 
is necessary for change in modifier frequency. More specifically, if we assume that the driver 
is costly and this cost is less than fully dominant (i.e. wdd ^ WDd ^ 1, and wdd < 1), an 
unlinked recombination modifier increases in frequency when it is underrepresented in drive 
haplotypes (i.e. LD < 0), as it avoids the fitness cost accrued by driving alleles. 

Analysis of the case of female-limited drive yields a similar conclusion. With sex-limited 
drive, the change in modifier frequencies in sperm and eggs equals 

A f a f f MtS ~ fM9 + Z /W (n i \ 

A /MdMI = A hi6Mll = 2 ^ ' 

a /a/9MI = a /j\/9MII = 2 ^ ' 

respectively, where 

z = (w D d - l)(LD<j/ d9 + LD 9 / dd ) + (wdd - w Dd ){LD<sf D9 + LD 9 f D6 ) (4) 



The first term in equations (13 .bp and (13. cp . (i.e. Jm6 — /a/9, and fuq — fnd, respectively) 



captures the role of syngamy in homogenizing allele frequencies across the sexes. The second 
term, z/W, is the change in modifier frequency due to linked selection on drive alleles. Since 
we assume that drive entails a less than dominant fitness cost (i.e. wdd ^ wod ^ 1, and 
wdd < 1), z will be positive when LD is negative. Therefore, recombination modifiers will 
increase in frequency when they are underrepresented on the driving haplotypes, like the case 
above. Additionally, since z/W has the same role in changing modifier frequency (A/m) in 
male and female gametes, this model does not generate a sex difference in modifier frequency. 

Since negative LD between recombination modifier and drive alleles is necessary for the 
adaptive evolution of the recombination rate, the salient question is - does our model generate 
negative LD? To address this question, we begin with a population in linkage equilibrium 
(LD = 0), and investigate the level of LD in gametes after selection, recombination, and 
drive. 
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Starting from LD = 0, the LD generated in a single generation with MI drive in both 
sexes equals 

LD M1 = Sr fDfdfAif m w Dd (l - 2a Ml )/W (5) 

The LD generated with Mil drive in both sexes equals LD'y^[ = — ^^MI' ^ e "^"^ generated 
between a female-limited driver and a recombination modifier after a single generation of 
selection, recombination, and drive is equivalent to the case described above replacing Sr 
and a by 5rg and ag. 

Note that with MI drive, negative LD is created between driving and recombination 
enhancing alleles (Equation ([5]) is less than zero when 5r and 5rg are positive). By contrast, 
with Mil drive, negative LD is created between driving and recombination suppressing alleles 
(when 5r and 5r<? are negative). Since negative LD between driver and modifier results in 
an increase in modifier frequency (e.g. Equations (3. a) and (4)), recombination enhancers 
increase in frequency by hitchhiking with high fitness non-driving haplotypes. Similarly, 
recombination suppressors are favored with Mil drive. 

We complement our single generation view of the creation of LD by making use o: 



the Quasi Linkage Equilibrium (QLE) method (e.g. 



Barton and Turellii . 



1991 



KlRKPATRICK et al. 



KlMURA 



1965 



Nagylaki 



1976 



20021 ) . QLE relies on the fact that when 



linkage is loose, and selection and drive are weak, LD approaches an equilibrium value on a 
much faster time-scale than the slow change in allele frequencies. To derive this QLE LD, 
we solve for the equilibrium level of LD* (i.e. the LD for which A LD = 0), while holding 
allele frequencies constant. To make this solution analytically tractable, we further assume 
that the strength of drive (a — 1/2 and a? — 1/2), selection (1 — wod and 1 — wdd)-, and 
that the degree of recombination modification (8r and 5r$), are all on order £, and we ignore 
terms of higher order than £ 2 . Doing this we find that 



LD* * 25rf d f D f m f M (l - 2a M i) 



(6.a) 



Comparing equation ( 16. aft to equation (jSJ) shows that much the equilibrium LD is generated 
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in a single generation. Recalling that by definition, aui > 1/2, it is clear that a recombination 
enhancer will be in negative LD with the drive allele, and will therefore be selectively favored. 
When considering Mil drive, the result is conceptually similar but of opposite sign - with Mil 
drive recombination suppressors, rather than enhancers generate negative LD with drivers. 

Employing QLE methods, we can also derive the equilibrium LD between recombination 
modifiers and female-limited drivers, with the same assumptions as above. In this case, the 
equilibrium LD in sperm is a function of the LD in eggs, LDJ = LDp/3. Combining this 
result and our finding that Jmq = fhid , and assuming weak selection, the equilibrium LD 
in eggs with MI drive is approximately 

T n * 3<5r(l - 2a)f Dd f Mm 

LDg ~ A-(l-2a)(l-r 9 )(f D -2f d6 ) ^ 

Again, the QLE results show that most of the LD in females generated in a single generation. 
Note that, although selection does not directly generate LD in males (above), the inheritance 
of haplotypes from females maintains a low equilibrium level of LD in males. 

In summary, when unlinked to a female meiotic driver, a female rec ombination modifier 



can spread by reducing the ability of the driver distort meiosis in its favor ([Thomson and Feldman 



1974 



Haig and Grafen 



19911 ). By decreasing the strength of drive, the recombination 
modifier becomes under-transmitted on the drive haplotype, as the driving allele drives less 
efficiently in the presence of the modifier. Therefore the modifier hitchhikes with the high 
fitness, non-driving allele; however, unlike traditional models of hitchhiking, the modifier es- 
sentially arranges a ride for itself by increasing the expected transmission of the nondriving 
allele. 

Furthermore, when drive is sex-limited, only recombination in the driving sex can act to 
decrease the efficiency of the single-locus drivers described above. Therefore, a modifier 
with equal and opposite effects of male and female recombination rates can spread, and so 



our results meet the criteria of 



LenormandI $2003) for the evolution of heterochiasmy. 



Predictions concerning the evolution of recombination modifiers unlinked to MI and Mil 
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A) Ml drive in both sexes B) Female limited, Ml drive. 



o 
d 




Generation Generation 



Figure 3: The convolution of an MI driver and an unlinked recombination en- 
hancer. The frequencies of MI drive alleles (fn, red), and unlinked recombination modifiers 
across generations. The correlation between alleles, ^7j=f=j^j ■> denoted by the 
red blue line, and its value is given on the right axis. Drive is complete and recessive lethal 
{wDd = 1) Wdd = 0). The initial recombination rate is 1/4, and each copy of M increases the 
probability of recombining by 0.05. Initial frequencies of drive and recombination modifier 
alleles equal f Do = 0.10 and /a/ = 0.01, respectively. [3]A) Drive in both sexes («mi = 1, 
r = 1/4, 5r = 0.05). 0p) Female-limited drive (a M i9 = 1, Sr 9 = 0.05). 

drivers are summarized by cases 1 and 2 in Table 1, respectively. Enhancers of the female 
recombination rate are favored when unlinked to female-limited MI drivers (case 1), while 
female recombination suppressors are favored when unlinked to female-limited Mil drivers 
(case 3). We display the dynamics of the case of the coevolution of a MI driver and an 
unlinked recombination modifier in Figure |3j 
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3.2 Single-locus drive - A linked recombination modifier. 

We now explore the fate of a recombination modifier in tight linkage with the drive locus 
(i.e. there is no recombination between drive and recombination modifier loci). Intuitively, 
a modifier linked to a driver faces two opposing pressures - the deleterious effect of drive on 
individual fitness, and the selfish effect of drive on allelic transmission. 

These two components are captured by the equation describing the change in modifier 
frequency when drive operates in both sexes: 

a r -LD(f d (l - w Dd ) + f D (w Dd - w DD )) LD w Dd (l - (r + 5r)) 

A /mMI = w + W (7 ' a) 

a p -LD(f d (l - w Dd ) + f D {wDd ~ Wdd)) , LDw Dd {r + 5r) 

A 4/MII = w + W ( } 

The change in frequency of a recombination modifier by individual selection is represented 

by the first term of equations (I7.al) and f!7.bj) and is equivalent to the unlinked case (equation 



(13. ap ). As in equation (13 ,ap . this term is positive when LD is negative. The later term in 
equations (17. al) and (I7.bj) represents the change in frequency of the modifier due to drive, 
which is positive when LD is positive. Thus, although individual selection favors linked 
recombination modifiers in negative LD with drivers, transmission distortion favors modifiers 
in positive LD with drivers. 

With female-limited drive, the intuition is similar; however, the recursive equations are 
more complex. The change in frequency of a recombination modifier in sperm, A/jy/d; is 
determined entirely by individual selection, and equals |(/jvfd — /mq + z/W), where z retains 
its value from equation (j4j). The change in the frequency of a recombination modifier in eggs 
(A/ M9 ) equals 



A / A /9MI = \Um9 fM6 + ^ + - (r 9 + 5r 9 ))) (7.c) 

A 4/9MII = \^M9 - Im6 + ^ + ^T^9 + M) (7-d) 
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where 

u = LDq + LD 6 + (f D9 - f D6 ){f M Q - /Md) (8) 

Therefore, with either sex-limited drive or drive in both sexes, both negative and positive 
LD between recombination modifier and driver can contribute to an increase in modifier fre- 
quency. Like the unlinked case, selection can generate linkage disequilibrium in a population 
in which no LD existed previously (see appendix). However, since a novel mutation must 
arise on one haplotype, the LD formed by the mutational history of tightly linked loci is 
more important than is the LD generated by selection. 

We display the population genetic dynamics of recombination modifiers linked to MI 
drivers in Figure HI Recombination enhancers increase in frequency when they arise on the 
non-driving background (Both sexes, Figure HJA.. Female-limited, Figure HB- See also Table 
1, case 1). By contrast, recombination suppressors spread when they arise on the driving 
background (Both sexes, Figure HP . Female-limited, Figure Hp. See also Table 1, case 2). 
The opposite result holds in the case of Mil drive (Table 1, cases 3 and 4) 

To provide a stronger intuition of the evolution of a recombination modifier linked to a 
drive locus, we investigate the special case of a recessive lethal driver which distorts meiosis 
in both males and females {wod = 1, wm = 0, and a = 1). In this case, Equations (I7.a[) and 
(I7.bp become 

A 4/MI = LD ( l -Id- (r + 5r))/W (9.a) 
A /mMII = LD ( r + Sr ~ fn)/W (9.b) 

respectively. Under the assumption of recessive lethality and complete drive, the equilibrium 
frequencies of MI and Mil drivers is straightforward and equals f^ m = 1 — r and fp MU = r, 
respectively (see appendix). Plugging these values into equations (19 .ap and (19. bl) we find 
that the change in frequency of recombination modifiers in tight linkage with MI and Mil 
drivers equals 

A f A iMl = -LD5r/W (9.c) 
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A/ mMII = LD6r/W (9.d) 

respectively. Equations (19. cp and (19. dp provide a straightforward characterization of the 
invasion of a rare recombination modifier. With MI drive and tight linkage (equation (19. eft ) a 
recombination enhancer will increase in frequency when it arises on the non-driving haplotype 
(LD < 0,Sr > 0, gives ^f M y[\ > 0), while a recombination suppressor will increase in 
frequency when it arises on the driving haplotype (LD > 0, 5r < 0, gives A/j^MI > 0) The 
opposite result holds for a recombination modifier tightly linked to an Mil driver (Equation 
( 19.cip ). We note that equations (19 .cp and (19. dp only hold for the first generation of selection. 

The analysis of female-limited drive is more complex, but ultimately yields a similar 
result. Here, we present our invasibility analysis in which we derive results under female- 
limited drive. In this case, female-specific recombination enhancers are favored when in 
repulsion phase with MI drivers or coupling phase with Mil drivers, while the opposite 
results hold for recombination suppressors (s ee 'Invasibility analysis ,' pages 18-19, 21 in our 



appendix). For this invasibility analysis (see 



Otto and Day 



20071 ). we write the recursion 



for haplotype frequencies in sperm and eggs after selection, recombination, and drive in 
matrix form. 

When eigenvalues of the Jacobian of this matrix (evaluated at viability-drive equilibrium 
setting the frequency of the recombination modifiers to zero) are greater than one, a rare 
recombination modifier increases in frequency. For both the cases of MI and Mil drive, 
there are two eigenvectors of interest. One eigenvector is greater than one when 5rg > 0, 
the other is greater than one when 5rg < suggesting that both recombination enhancers 
and suppressors can increase in frequency). Because the LD generated in these cases is 
large (in fact, at equilibrium there are only two haplotypes), we do not e mploy the QLE 



approach, which works best under loose linkage 



1965 



KlRKPATRICK et al. 



2002 



Nagylaki 



1976|). 



Barton and Turelli 



1991 



KlMURA 



Unlike the case of unlinked recombination modifiers, modifiers linked to drive loci do 
not generally approach fixation. Assuming that the modifier is favored, it rapidly goes to 
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fixation on the background onto which it mutated; however, this haplotype now moves to its 
new equilibrium frequency determined by the new recombination rate. So long as this new 
equilibrium is greater than zero and less than one, recombination modifiers will be stably 
polymorphic (see Figure H]). 

In sum, the evolution of recombination modifiers linked to drivers yields a rich and diverse 
set of predictions. When tightly linked to MI drivers, recombination enhancers in repulsion 
phase (Table 1, case 1), and recombination suppressors in coupling phase are favored (Table 
1, case 2). When tightly linked to Mil drivers, recombination suppressors in repulsion phase 
(Table 1, case 3), and recombination enhancers in coupling phase are favored (Table 1, case 
4). Like the unlinked model (above), the fate of a female recombination modifier linked to a 
female-limited driver is independent of its influence on the male recombination rate. Thus, 
a modifier with equal but opposite effect on male and female recombination rates (i.e. no net 
effect) can spread, facilitating the evolution of heterochiasmy. 



4 Two-locus drive systems 



We now turn our attention to the more complex case of two-locus, MI drive. In this model the 
strength of drive by a centromeric variant, C, depends on the genotype at the drive-enhancer 
locus, D, which is on the same chromosome as the centromeric driver. Specifically, in Cc 
heterozygotes, meiosis is fair in a genetic background of d homozygotes, but C is represented 
in ct\ and cti of gametes from Dd/Cc and DD/Cc individuals, respectively (where cti ^ ct\ ^ 
0.5 and cti > 0.5). Although it is possible that the drive enhancer will incur an individual 
fitness cost, we focus on the case in which the drive enhancer is neutral, but the driving 
centromere is costly. Imposing a fitness cost to drive enhancer adds subtle quantitative 



differences to the results, and this mode 



in the description o : 



Haig and Grafen 



the SD system (IHartl 



where b oth loci involve costs has been well explorec 



1975 



Charleswqrth and Hartl 



1978 



19911 ). Since the genetic identity of a centromere seems unlikely to 
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Figure 4: The evolution of drivers and recombination modifiers in tight linkage 

The frequencies of MI drivers (/d, red), and linked recombination modifiers (Jm, blue) across 

generations. The correlation between alleles is denoted by the red blue line, and its value 

is given on the right axis. Initial frequencies of driver and recombination modifier alleles 

are fjj = 0.10 and /m = 0.01, respectively. Drive is complete and recessive lethal. Initial 

recombination rate equals 1/4. [4]A.) Drivers and recombination enhancement in both sexes 

(Sr = 0.05), M arises on a d chromosome. |4jEJ) Female-limited driver and recombination 

enhancement. (5rg = 0.05), M arises on a d chromosome. |4]C) Drivers and recombination 

suppression in both sexes. (8r = —0.05), M arises on a D chromosome. HP) Female-limited 

driver and recombination suppression. (5rg = —0.05), M arises on a D chromosome. 
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influence Mil drive, we do not pursue a two-locus model of Mil drive. 

With two-locus, MI drive, a recombination enhancer can incr ease in frequency, and ulti 



mate l y approach fixation, as in t he case of single- locus MI drive ([Thomson and Feldmani . 



1974 



Haig and Grafen 



19911 ). With drive in both sexes, the change in frequency of a 



recombination modifier is 

Af M = ^^(f c (l - w Cc ) + fc(w Cc ~ wco)) (2) 

where LD MC = juc — fufc is the linkage disequilibrium between centromere and recombi- 
nation modifier. 

As in the case of recombination modifiers unlinked to single-locus MI drivers, the re- 
combination enhancer spreads by becoming underrepresented in the low fitness, centromeric 
driving (C) genetic background. Recombination enhancers generate this LD by decreas- 
ing the expected transmission of the drive enhancer allele, which allows the recombination 
enhancer to escape from the driving haplotype (Figure HJA.). 

When drive is female-limited, alleles that increase the recombination rate between drive 
enhancer and centromeric loci in either sex are favored. However, female-limited recom- 
bination enhancers spread much more quickly and have more negative LD with driving 
centromeres than do male-limited recombination modifiers. For example, in Figure HJ3 it 
only takes approximately 11,400 generations for a female-limited recombination modifier to 
rise from a frequency of 0.1% to 95%, but it takes more than an order of magnitude longer for 
a male-limited recombination enhancer to reach this frequency (Note the order of magnitude 
difference on the x-axis in [4)3 and HP). 

Because female recombination enhancers are more strongly favored than male recombi- 
nation enhancers in this system, alleles that increase female recombination can increase in 
frequency, even if they drastically reduce male recombin ation rates (Figure HP ). Therefore, 



our two-locus model also passes the 'no- net-effect test' (ILenormand 
the evolution of heterochiasmy. 



20031 ) - facilitating 
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The intuition behind this numerical result is as follows. Drive in females generates a pos- 
itive association between centromeric drivers (C) and linked drive enhancers (D). With this 
association, D increase in frequency during female (but not male) meiosis. Recombination in 
females decreases both the expected transmission of drive-enhancers, and the co-transmission 
of D and C alleles. Since gametes from mothers with higher recombination rates have fewer 
D alleles than expected, drive in their daughters is less efficient. Ultimately, the granddaugh- 
ters of females with hig her recombina tion rates suffer less from the deleterious fitness effects 



of drives systems (e.g. 



Crow 



19911 ). 

As male recombination does not directly change the transmission of a female meiotic drive 
enhancer, the selection on a male recombination modifier is weak. Nonetheless, elevated male 
recombination rates do break down the association between D and C alleles, which ultimately 
allows male recombination enhancers to escape from centromeric drivers, providing a minor 
boost to male-limited recombination enhancers. However, this effect is meagre compared to 
the effect of female recombination modification (Figures HB-D). 
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Figure 5: The co-evolution recombination modifiers and a two-locus drive sys- 
tem. The frequencies of drive enhancer, centromeric driver, and recombination modifier 
alleles are displayed in red, green, and blue, respectively. The correlation between recom- 
bination modifier and centromeric driver alleles is denoted by the green blue line, and its 
value is given on the right axis of Figure 5B (The scale is maintained in Figure 5C). Ini- 
tial frequencies drive and recombination modifier alleles equal fo = 0.10 and /m = 0.01, 
respectively. The driving centromeric allele completely distort meiosis in DD and Dd geno- 
types (e.g. oji = 0J2 = 1), and is a recessive lethal. Neither drive enhancer nor recombination 
modifier directly influence individual fitness. The initial recombination rate and allele fre- 
quencies are: r = 0.10, and fn(o) = fc(o) = /m(o) = 0.001, respectively. A) Recombination 
modification in both sexes (5rg = 0.025, 5ra = 0.025). B) Female-limited recombination 
modification (dry = 0.025, dr^ = 0). C) Male-limited recombination modification (5r$ = 0, 
5rg = 0.025). D) The modifier has distinct influence on male and female recombination 
rates (Sr^, and 5rg, respectively). Green indicates an increase in modifier frequency, pur- 
ple indicates a decrease. Labels above diagonal lines describe the relative change in allele 
frequencies over 250 generations. 
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5 Discussion 



Meiosis and recombination are deeply conserved and highly orchestrated processes, in which 



slight errors can have severe fitness consequences. Nonet heless, many o : 



ponents of meiosis and reco mbination evolve rapidly (e.g. 



Malik et al. 



the functional com- 



2002 



Anderson et al. 



2009 



Myers et al 



2010j). One explanation for this rapid evolution is that meiosis and ga- 



metogenesis offer a number of opportunities for genomic conflict within an individual, gener- 
ating a pattern of a ntagonistic coevolution betw een selfish gametic drivers and suppressors 



of meiotic drive (see 



Burt and Trivers 



2006. for a broad overview). 



Since the progression of meiosis and gametogenesis is highly sex-specific (IMorelli and Cohen 



20051 1. we expect that forms of conflict and conflict mediation will also be sex-specific. The 
asymmetry in meiotic division during oogenesis presents an opportunity for comp e tition 



Zwick et al. 



1999 



alleles for representation in the egg ( 


Sandler and N 


OVITSKI 


1957 


Pardo-Manuel De Villena and Sapienza. 


2001a 


b 


Malik. 


2009) 



Because recombination determines the ability of female drivers to distort meiosis at MI and 
Mil, female recombination rate influences the ability of a driver to distort meiosis. We have 
shown that female meiotic drive can favor changes in the female recombination rate - female 
recombination modifiers are selected to enhance or suppress drive (see Table 1) with changes 
in male rate having little or no effect on the efficacy of female drive. 

Our understanding of the frequency, severity, and operation of female meiotic drive sys- 
tems is still in its infancy. Therefore, we have not based our population genetic analysis on 
explicit mechanistic details of female drive. However, our models can be related to biolog- 
ically plausible mechanisms. The model of a single-locus MI driver could correspond to an 
epigenetic modification of a centromere in cis, or a structural stretch of DNA that influences 
the orientation of the centromere in such a way as to increase its probability of inclusion 
into the secondary oocyte. In our two-locus MI drive system, the centromeric locus could 
correspond to a centromeric satellite that increases its probability of inclusion in the primary 
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oocyte through an interaction with the spindle, while alternative alleles at the drive enhancer 
locus could represent centromeric protei ns which interact with the centromeric machinery to 



enhance or su ppress the effect of drive (IMalik et al. 



Malik 



2002 



Malik and Henikoff 



2009 



20091 ). Our Mil mod el roughly corresponds to neocen tromeric drive systems, such 



as the Ab-10 locus in maize (IRhqades and Dempsey 



the orientation of meiotic chromosomes ( 



NOVITSKI 



1951 



sec 



1966[). or telomeres that infl uence 



Anderson et al. 



(J2008|) for a 



discussion). We reiterate that our results depend solely on the ability of recombination to 
modify drive systems, rather than on any specific drive mechanism. 

We have shown that female meiotic drive systems create a selective environment that 
favors the evolutionary modification of the female recombination rate. However, selection 
on female recombination rates does not necessarily lead to heterochiasmy. For example, in 
the extreme case where if alleles that modify the recombination rate in females have equiv- 
alent effects on the male rate, selection on the female rate will not generate heterochiasmy. 
However, above we show that our model favors modification of the female recombination 
rate even if modifiers have o pposing effects the o n male rate - the standard for models of the 



evolution of heterochiasmy ( ILenormand 



20031 ). 



While we still know little about recombination modifiers, current evidence suggests that 
the genetic control of the global recombination rates differs by sex, and therefore selection on 
the female rate is likely to generate heterochiasmy. Three distinct lines of evidence support 
this tentati ve conclusion. First, the glob al control of meiosis and recombination is sexually 



dimorphic (IMoRELLl and Coheni . 



20051 ). Second, the few naturally occurring alleles known 



to modify th e tota 



( Kong et al 



2004 . 



genet i c map length in humans d o so in a highly sex-specific manner 



2008 



Fledel-Alon et al 



201 ll ). Third, although there is additive 



genetic variation for the ma p length in both sexes, no heritable intersexual correlation in 



map length has been found (IFledel-Alon et al. 



20 111 ). 



The predictions of our model are sensitive to biological details such as the linkage asso- 
ciation between drivers and recombination modifiers, as well as the timing of drive (MI vs 
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Mil). This makes it hard to generate concrete predictions about whether female drive will 
select for higher or lower female recombination rates. Indeed, the fact that female rates are 
not always higher than male rates suggests that the direction of selection may not be con- 
stant. One concrete prediction is that since the centromere holds a special place in female 
meiotic drive (in both MI and Mil systems) a constant influx of new female drivers will 
systematically select particularly for heterochiasmy in the regions surrounding centromeres. 

The observation in many taxa of higher female recombination rates, especially near the 
centromere is consistent with two different biological scenarios elaborated in our models. 
First, elevated female recombination rates, especially near centromeres, may represent the 
action of unlinked suppressors to prevent the spread of MI female drivers. Alternatively, 
this pattern could be explained by the spread of recombination enhancers linked to Mil 
drivers, which increase in frequency because recombination enhancers facilitate Mil female 
drive. Empirical progress in elucidating the genetic basis of local variation in sex-specific 
recombination rates and female meiotic drive across the genome will shed light on which (if 
any) of these two models explain this pattern. 

In contrast to global modifiers of the recombination rate, we know a little more about 
local modifiers of recombination rate, though our picture is still far from complete. One 
broad class of local suppressors are chromosomal in versions, which seem to be a common 



KlRKPATRICK 



2010 



for discussion). 



response to selection for reduced recombination (see 
Inversions are a priori expected to locally suppress recombination similarly in both sexes 
and this heterozygous effect will be removed when the inversion is eventually lost or fixed 
within the population. Therefore, given our current understanding of local recombination 
modification, we think it is unlikely that selection for linked suppressors of recombination 
will strongly contribute to heterochiasmy. 

Although we focused on female-limited drive, there are well documented cases of male- 
limited transmission distortion involving multilocus gen e complexes (Note that no ne of these 



are true meiotic drive, relying instead on sperm death) ( jWu and Hammer 
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199 lh . A model 



of the coevolution of a two-locus male drive system and a recombination modifier will yield 
results similar to the case of two-locus MI female drive described above (with slight differences 
due to the fitness of recombinants in male systems). That is, with male-limited drive systems, 
we expect that recombination modifiers in coupling phase would benefit from reducing the 
male recombination rate, while unlinked modifiers would benefit from increasing the male 

2~ 



recombination rate (see also discussion by 



Lenormand and Dutheii 



20051 ). Evidence for 



the former is bountiful, a s most known male transm ission ratio distorters are tied together by 



complex inversions (e.g. 



Presgraves et al. 



20091 ). However, unlike female meiotic drivers, 



male distortion systems can arise at any chromosomal location irrespective of the distance 
to the centromere. Therefore, even if male distortion systems do select for heterochiasmy, 
a constant influx of new male systems will not systematically select for sex differences on a 
chromosomal-scale. 

More generally, our models may explain other broad patterns associated with recom- 
bination. One outstanding pattern is the observation that variation in the number of 
chromosome arms is a better predictor of recombination rate variation among mammalian 
species than is variation in the number of chromosome s, or the physical size of the genome 



(IPardo-Manuel De Villena and Sapienza 



200ld ). This result is somewhat surprising 



because only one recombination event per chromosome is required for p roper segregation, 



arms of metacentric chromosomes are often found to be lacking a crossover (IFledel-Alon et al 



2000 



2009 ) , and the centromere seems to offer no barrier to interf erence in many systems ( 



Fledel-Alon et al. 



2009 



Broman and Weber 



Demarest et al. 



20111 ) . The meiotic drive theory can ex- 



plain these observations by proposing that modifiers of the recombination rate are selected 
to increase recombination events between centromeres and potential drivers, as both sides 
of a chromosome present an opportunity for drivers to exploit meiosis. 

An other broad pattern is the ob s ervat ion that heterochiasmy is reduced in selfing plants, 



which 



Lenormand and DutheilI (120051 ) saw as support for their hypothesis of a role of 



haploid selection on male gametes in driving heterogamy. We note that this observation is 
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potentially consistent with the female meiotic drive hypothesis - since most selfing plants are 
largely homozygous, there is little opportunity for drive. The alteration of generations in 
plants provides exciting opportunities for future research on the evolution of heterochiasmy. 



The meiot ic drive theory could also 



nation rates (ICOOP and Przeworski 



explain the observation of rapic 



2007 



Dumont and Payseur 



changes in r ecombi 



2008 



2011( 1. as the 



recombination map will constantly be evolving as recombination enhancers and suppressors 
respond to new drive systems across the genome. A greater knowledge of medium-scale pat- 
terns of turnover in male and female rates would help clarify the plausibility of meiotic drive 
in explaining this pattern. For example, the meiotic drive theory would be strongly sup- 
ported if regions close to centromeres show parti cularly high rates o f female recombination 



19961). 



evolution, as is observed in the Drosophila clade ( ITrue et al. 

Further tests of these predictions require the ability to identify both meiotic drivers and 
modifiers of the recombination rate. Currently, our knowledge of the distribution of meiotic 
drivers and their fitness effects is very incomplete, and is likely biased towards the overrepre- 
sentation of strong drivers with extreme fitness effects. However there is mounting ev i dence 



supporting the e xistence of subtle tran s mission ration distorter s (e.g 



Zollner et al. 



2004 : 



Aparicio et al. 



2010 



Axelsson et al. 



2010 



Reed et al. 



2005 



however, the mech- 



anism of distortion in these cases is often unknown). Similarly, although there is ample 
evidence that the recombination rate, as well as the strength and direction of heterochi- 
asmy varies across species, few allelic variants that influence sex-specific recombination rates 
have been identified. Although we currently know very little about the frequency of female 
drivers, or the genetic control of sex-differences in the recombination rate across many taxa. 
Fortunately, as technological advances make the sequencing of offspring (or gametes) from 
many meioses more affordable, identification of alleles governing sex-specific transmission 
and recombination rates will become much easier. 

More broadly, the ideas presented here are part of a larger body of the theoretical 
work that highlights the potential diverse role of genomic conflict in shaping the evolu- 
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Haig and Grafen. 


1991: 


Haig 


2010 




Hurst et al. 


. 1996 


2009; 


Anderson et al. 


2009 


Burt and Trivers 


■ * 


!006 





Sandler and Novitski. 



Zwick et al. 



1999 



1957 



Malik and Henikoff . 



, Table 12.3). For example, conflicts 



between cis and trans determinants of hotspot localization due to bias gene conversion have 



been put forward to explain the rapid evolution of mammalian fine-sea 



e recombination rates 



and their determinants (such as the hotspot bi nding protein Pr d m9 (IBaudat et al 



Myers et al. 



2007 



2010 



Parvanov et al. 



Ubeda and Wilkins 



2010|)) (Iboulton et al. 



1997 



2010 



Coop and Myers 



20111 ) . While much of the classic work on the evolution of re- 
combination and meiosis has focused on the benefits of creating adaptive gene combinations, 
purging deleterious recessive alleles from an adaptive haplotype, or bringing together two 



beneficial mutations onto one haplotype (e.g. 



1996 



Otto and Lenormand 



2002 



Es 



Barton 



hel and Feldman 



1970 



Feldman et al. 



20091 ). it seems equally plausible that the 



short-term evolution of recombination rates may be in response to conflicts created during 
meiosis. 
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